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. Abstract. I review the arguments motivating models for massive star 

. ' formation via stellar collisions. I then describe how the standard accre- 

. tion scenario, involving the collapse of a quasi-hydrostatic gas core, can 

produce high-mass stars in the pressurized regions of forming star clus- 
ters. I argue that the observational evidence, particularly in the Orion 
hot core, favors the standard accretion paradigm. 

m 1 

! 1. Introduction 

o 

Two basic models for how a massive star accumulates its mass are debated. The 
conventional case, which I refer to as the standard accretion model, involves the 
^S) \ inside-out collapse of a gravitationally bound, pre-stellar gas core from approxi- 

mate hydrostatic equilibrium. A small amount of angular momentum creates an 
accretion disk, via which mass reaches the star in an ordered manner. Low-mass 
Q H i stars appear to form in this way (Shu, Adams, & Lizano 1987). 



o 



The collisional model involves the coalescence of smaller stars or protostars 
in dense stellar clusters (Bonnell, Bate, & Zinnecker 1998). These lower-mass 
stars may first form via standard accretion from small gas cores that have frag- 
mented from the protocluster gas clump. They may also build up their mass via 
Bondi-Hoyle accretion from the clump. 

2. The Case Against Standard Accretion 

Some of these points are discussed in the review by Stahler, Palla, &; Ho (2000). 

Long Formation Time: Assuming 100% efficiency, the accretion rate 
of a star forming from the collapse of a critically unstable isothermal sphere is 
= 0.975c4/G = 4.4 x 10~ 6 (T/20K) 3 / 2 M yr" 1 (Shu 1977), so the formation 
time is t*f = m*f/rh» = 6.9 x 10 6 (T/20 K)~ 3 / 2 (m*//30 M & ) yr. This timescale 
should be less than the main sequence lifetime (several Myr for massive stars) 
and the cluster formation time (~ 1 Myr for the Orion Nebula Cluster (ONC); 
Palla & Stahler 1999), requiring large accretion rates and temperatures (e.g. 
m* > 10~ 4 M y^ 1 and T > 160 K for m*/ = 100 M and U f < 1 Myr). Such 
temperatures are observed in massive star-forming regions, but are unlikely to 
occur until at least one massive star is present. Higher accretion rates result if a 
core has some nonthermal pressure support (Stahler, Shu, & Taam 1980). This 
does appear to be the case for massive cores since their line-widths correspond 
to supersonic velocities. 
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Radiation Pressure: The short Kelvin-Helmholtz times of massive stars 
allow them to reach the main sequence while they are still accreting (e.g. Palla 
& Stahler 1992). These protostars have high luminosity, L, and the radiation 
pressure may disrupt infall by acting on dust grains, which survive until tem- 
peratures of ~ 2300 K and are well-coupled to the dense gas. A necessary 
condition for infall is that the ram pressure of inward motion of the gas ex- 
ceed the outward radiation pressure from the direct radiation field of the star at 
the dust destruction front, r d ~ 13(L/10 5 L o ) 1 / 2 (T d /2300K)- 2 (Q (l /0.1)- 1 / 2 AU, 
where Q a is the effective absorption efficiency of dust grains at the sublima- 
tion temperature. This implies pv 2 > L/(4 / Kr 2 i c), where p and v are the den- 
sity and infall velocity of the gas at r^. In the spherical case we then have 
m > L/(cv). For zero age main sequence (ZAMS) stars of mass 30, 60, 120 M , 
L ~ 1.2, 5.1, 18 x 10 5 L© (Schaller et al. 1992), and assuming free-fall conditions 
at rd we have m ^ 0.4,1.7,6.2 x 10~ 4 M© yr _1 (see also Wolfire & Cassinelli 
1987). Nakano (1989) and Jijina & Adams (1996) have argued that these condi- 
tions are relaxed for accretion via a disk, since once gas reaches the disk (typically 
at radii S> r^), it is shielded from much of the stellar radiation. However, the ram 
pressure criterion should still apply in the disk at r^. Consideration of this case 
(Tan & McKee 2002, in prep.) suggests that accretion rates ~ 10~ 3 M yr _1 
are required to form the most massive stars. Radiation-hydrodynamic simula- 
tions (e.g. Yorke & Sonnhalter 2002) will ultimately provide the most accurate 
answer to this question. 

Crowding: It appears that massive stars form almost exclusively in clusters 
containing many more low-mass stars, which dominate the total stellar mass. 
The central stellar densities can be very high — Hillenbrand &z Hartmann (1998) 
find n* ~ 1.7 x 10 4 pc -3 in the ONC corresponding to separations of about 0.04 
pc. The hydrostatic cores invoked in standard accretion would have to form and 
survive in such a crowded environment. 

Small Jeans Mass: The Jeans mass (or more appropriately the Bonnor- 
Ebert mass, which is the maximum mass of a stable isothermal sphere) becomes 
very small in regions of high pressure, such as high-mass star-forming clumps 
where P ~ 10 8 - 10 9 K cnr 3 (McKee & Tan 2002b [MT02b]). The Bonnor- 
Ebert mass is M BE = 0.046(T/20 K) 2 (P/10 9 Kcm- 3 )' 1 /^. If gas forms 
bound structures only on this scale, then to form a massive star, most of the 
mass would have to be accumulated via Bondi-Hoyle accretion or mergers. 

3. The Case for Collisions and Competitive Accretion 

A model of massive star formation by collisions of lower-mass stars (Bonnell et al. 
1998), whose masses may be augmented by competitive Bondi-Hoyle accretion 
(Bonnell et al. 2001a), can automatically circumvent three of the four difficulties 
faced by standard accretion. Firstly, the coalescence of stars, or any optically 
thick mass unit, is not prevented by radiation pressure. However, Bondi-Hoyle 
accretion is suppressed, i.e. for m* ;> 10 M©. Rates of Bondi-Hoyle accretion 
are enhanced in dense regions and collisions are favored in crowded regions — just 
the locations where massive stars are observed to form. Bonnell k, Davies (1998) 
have shown that there has not been enough time for the Trapezium stars in the 
ONC to have reached their central location by dynamical relaxation. Either 
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they formed in the center (half of the massive stars forming within 30% of the 
half-mass radius) or they formed with much smaller velocity dispersions than 
the lower mass stars. Bondi-Hoyle accretion and mergers can potentially build 
up sub-solar Bonnor-Ebert masses to the massive star regime. A simulation 
of competitive accretion (Bonnell et al. 2001b) has reproduced the observed 
stellar IMF from an initial cluster of a 1000 identical 0.1 Mq stars embedded in 
a cold gas clump that makes up 90% of the mass. Note, however, that in some 
star-forming regions pre-stellar gas cores are observed to have a mass spectrum 
similar to that of stars (up to at least 5 Mq) (e.g. Motte et al. 2001). 

It has been proposed that run-away OB stars, which are ejected from clus- 
ters with velocities ~ 50 — 150 km s , and the small primary to secondary mass 
ratios of massive star binaries, may be the by-products of massive star formation 
by collisions (Stahler et al. 2000). Both phenomena are expected to result from 
close three body interactions of massive stars — the least massive star tends to be 
ejected leaving the remaining stars in a tightened binary. Such encounters should 
be common if massive star mergers are occurring. An alternative mechanism for 
producing run-away OB stars (from clusters that are at least several Myr old) is 
a supernova in a double massive star binary system. Hoogerwerf, de Bruijne, & 
de Zeeuw (2001) have identified the ejected stars from single examples of both 
of these types of event. 

It is clear that collisions and strong dynamical interactions can occur in the 
crowded environments of star clusters, but the real question is whether the rate 
is great enough to be relevant for the star formation process. The collisional 
timescale in a cluster of equal mass stars in the strong gravitational focusing 
limit is (Binney & Tremaine 1987) 

— 1 — (^) (»=f0 (^) (*) * « 

where a is the ID velocity dispersion, and r* and m, the stellar radius and 
mass. In the above equation we have normalized the variables to values typical 
of the central region of the ONC and allowed for a generous increase in the 
stellar radius due to pre-main sequence activity. For collisions to be important 
we require the collisional timescale to be ^ 10 6 yr, which requires a substantial 
change in one or more of the above parameters. 

Stahler et al. (2000) pointed out that the collisional cross section is sub- 
stantially increased for protostars that are still embedded in dense gas cores or 
have massive disks (see also McDonald & Clarke 1995). However, it is also pos- 
sible that such encounters separate the stars from their gas without leading to a 
collision (Price & Podsiadlowski 1995), which would reduce the efficiency of this 
mechanism. Zwart et al. (2001) have found that mass segregation and binarity 
can boost the collision rate. It may also be increased if stars are swollen by the 
energy release from a recent collision. Nevertheless, even for an optimistic colli- 
sion radius of 1AU = 215R0 (for each star) and for stellar masses m* = 10M@ 
equation ([[]) shows that densities of order 10 6 pc -3 are needed for t co n ~ 10 6 yr. 
More realistic cross sections suggest the densities need to be ^ 10 8 pc~ 3 . 

Bonnell et al. (1998) have presented a model in which extreme densities 
can result in a cluster of lower-mass stars that are accreting from their proto- 
cluster gas cloud. The velocity dispersion of the stars decreases as they accrete, 
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lowering the total energy of the system and causing the cluster to contract. In 
a subsequent version of this model (Bonnell 2002) the initial gas cloud, which 
dominates the total mass, is in free-fall collapse. The cloud is seeded with var- 
ious spatial distributions of low-mass stars. Gas densities increase rapidly at 
the center and low-mass stars that happen to be here grow quickly with large 
rates of Bondi-Hoyle accretion (Bonnell et al 2001a,b). After one free-fall time 
the central stellar density can increase from its initial value by factors of about 
10 5 . However, this dramatic increase appears to depend on the assumption of 
free-fall collapse of the entire gas cloud. It is argued that if the initial conditions 
for this simulation are taken to be the present-day densities of the ONC, then 
collisional massive star formation is likely to occur. 

Finally, it has been suggested that filamentary structure in the collapsing 
gas could provide a means of weakly collimating any outflows that are produced 
in the star formation process (Bonnell 2002). It would be difficult for the colli- 
sional model to produce coherent and symmetric outflows. 

Two key observational predictions of the collisional model are the sporadic 
release of large amounts of energy from stellar collisions and the brief existence 
of dense (^ 10 6 — 10 8 pc -3 ) stellar clusters around each forming massive star. 
Neither of these have been observed, although they may be obscured by the high 
extinction (Ay ^ 200) typical of massive star-forming regions. 

4. Turbulent Core Accretion 

In spite of the difficulties facing the standard accretion model (§2), the presence 
of coherent and massive gas cores in high-mass star-forming regions (Garay &: 
Lizano 1999; Kurtz et al. 2000) suggests that a model for massive star formation 
should start by considering the collapse of such structures (McLaughlin & Pu- 
dritz 1997; Osorio et al. 1999; McKee & Tan 2002a [MT02a] ; MT02b). Massive 
stars do not form in isolation and in the model of MT02a,b we assume the cores, 
which may form individual or binary stars, are part of a self-similar hierarchy of 
structure. We term the gas cloud that forms a star cluster a clump. The density 
distributions of clumps and cores, averaging over internal dumpiness, are consis- 
tent with approximately spherically symmetric power law density profiles, with 
p o<- r -k P an j £ ^ 15 (Evans et al. 2002). Such structures can result from the 
hydrostatic equilibrium of gas with an equation of state P oc p 7p , with 7 P = 2/3. 
Since the gas is observed to be cooler than the virial temperature, it must be 
supported by nonthermal forms of pressure. The signal speed, c = (P/p) 1 ^ 2 , is 
therefore supersonic and so the cores and the clumps should be turbulent and 
clumpy, as is observed. In the fiducial case of k p = 3/2, c oc r 1 / 4 and P oc r~ kp , 
with kp = 1. We define the surface of a core to be where the pressure has de- 
creased to the ambient pressure in the clump. This is set by self-gravity and is 
P c i ~ 0.88GS^ = 8.5 x 10 8 S^K cm -3 in the central regions where massive stars 
typically form. In this expression S c i = M c \/(itR 2 ^) is the mean surface density 
of the clump, with characteristic value ~ 1.0 g cm -2 (Plume et al. 1997). 

Since the core mass inside a distance r from its center is M(< r) = kpc 2 r/G, 
then for higher ambient pressures (larger values of c 2 at the core surface), the 
equilibrium state for a given core mass becomes smaller and denser. The ef- 
ficiency (e cor c = m*f /M corc ) of star formation from a core is typically quite 
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high for low-mass stars (~ 50%-Matzner & McKee 2000), if set by magneto- 
centrifugally driven bipolar outflows such as disk-winds or X-winds (Kdnigl & 
Pudritz 2000; Shu et al. 2000) that launch a fraction f w ~ 0.1 - 0.3 of the 
accreted mass with force p w = m w v w = f w m*v w = 4> w rh^VK- Here vk is the Ke- 
plerian velocity at the equatorial stellar radius and 4> w = 0.6 is a fiducial value. 
These winds sweep up and eject gas from polar regions of the core. Similarly 
high efficiencies hold for massive stars for the same values of f w and cj> w (Tan & 
McKee 2002, in prep.). For e core = 0.5 we have (MT02b) 

^"MsSfey'V'V' <2) 

Recall the mean stellar separation in the center of the ONC of about 0.04 pc. 
Comparison with the expected size of a massive core shows that crowding is not 
necessarily a problem, particularly since during a typical stage of formation half 
the stars observed today were not yet present, in the central part of the ONC 
there has been time for some relaxation and density increase since formation, 
and most stars are low-mass and have little dynamical impact on a massive core. 

When a core does collapse it must be from a configuration at least as dense 
as implied by equation (|2[). The high pressures of massive star- forming clumps 
require dense cores, which have short free-fall times. Inside-out collapse proceeds 
as in the model of Shu (1977), but now, since c oc r 1 / 4 , the accretion rate grows: 

where m* is the instantaneous protostellar mass. It will fluctuate because of 
core dumpiness. These mean accretion rates are large enough to overcome 
the radiation pressure of massive stars (Tan & McKee 2002, in prep.). The 
corresponding formation timescales, 

are short compared with the age constraints of the ONC (Palla & Stahler 1999). 

In these models the mass of a star is determined by the initial core mass 
and the efficiency of star formation, e core - Since we find e core to be relatively 
constant with stellar mass, the shape of the stellar IMF should reflect that of 
the cores. This is consistent with observations (e.g. Motte et al. 2001). Such a 
mass function may result from the coagulation of cores or from fluctuations in 
the turbulent velocity field that produce unstable condensations, but we have 
not attempted to predict it. We do not expect much hierarchical fragmentation 
of collapsing cores if the collapse happens quite quickly after formation (i.e. 
within a few dynamical times) and if most turbulent fluctuations do not produce 
unstable cores. The central concentration of the cores (p oc r~ 15 ) may suppress 
fragmentation (Bodenheimer et al. 2000). Once a protostar has formed, its tidal 
field will also stabilize the collapsing core against fragmentation. 

We have calculated the evolution of protostellar radius and luminosity (bolo- 
metric, outflow-mechanical, and ionizing) for the accretion rates predicted from 
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Figure 1. Models of protostellar evolution: bolometric luminosity 
(top panel) , outflow force (middle panel) , and ionizing luminosity (bot- 
tom panel) versus protostellar mass. Solid lines show model predictions 
for m*f = 7.5, 15, 30, 60, 120 M stars forming in a S c i = lg cm -2 
clump. Dashed and long-dashed lines show a 30 M star forming 
from S c i = 0.32, 3.2g cm~ 2 clumps respectively. The dotted line shows 
the ZAMS bolometric (Schaller et al. 1992) and ionizing luminosities 
(Schaerer & de Koter 1997; W. Vacca, private comm.). Constraints on 
the properties of several nearby protostars are derived from the lumi- 
nosities of their hot molecular cores (see Osorio, Lizano, & D'Alessio 
1999 for additional modeling and observational references). 
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the turbulent core collapse model. Before reaching the ZAMS, the radius is de- 
termined by deuterium core and shell burning (e.g. Palla & Stahler 1992). Our 
model is based on that of Nakano et al. (1995), but extended to include shell 
burning with a simple prescription that matches the results of Palla & Stahler. 
With the protostellar radius and accretion rate we calculate the accretion lumi- 
nosity, which adds to the internal luminosity to give the total. The outflow force, 
p w , is calculated assuming (fi w — 0.6, as described above. The rate of emission 
of H ionizing photons, S, is calculated with contributions from the star and the 
accretion shock and utilizing models of Schaerer & de Koter (1997) for solar 
metallicity. The models are shown in Fig. 1. From the observed bolometric 
luminosities of hot molecular cores, we constrain masses and accretion rates of 
embedded protostars, assuming S c i = lg cm -2 . This also constrains p w and S. 

If the mechanism of magneto-centrifugal outflow generation remains quali- 
tatively similar to that operating in lower-mass stars, where theoretical models 
(Shu et al. 2000; Konigl & Pudritz 2000) predict f w = m M /m* ~ 0.1 - 0.3 
with most of the outflow originating from the inner (r <J 10r*) accretion disk, 
then the large accretion rates predicted by equation (||) imply a high density of 
the outflowing gas rising from the disk. For fiducial parameters, we find that 
the ionizing fluxes of protostars at least as massive as ~ 60M@ are confined 
by their outflows in directions approximately in the plane of the accretion disk. 
Thus we do not expect disk photoevaporation (Hollenbach et al. 1994) to be an 
important process during the accretion phase of most massive stars. 

Ionizing photons can more readily escape along the polar directions of the 
flow, where MHD-outflow models predict the existence of "dead-zone" cavities 
that the outflow does not penetrate. These regions should be occupied by a more 
diffuse conventional stellar wind. The angular extent of the cavity as viewed from 
the star decreases with distance as the outflow is collimated by the hoop stresses 
of the toroidal magnetic field component. Outside the cavity, the density of the 
outflow may be derived given the angular distribution of the outflow force and 
assuming a constant terminal velocity. For the force distribution of Matzner & 
McKee (1999), which includes the small angle softening parameter 0q, we have 

4.34 x 10 5 f f w 10AU \ 2 0.6 ln(200) m« / 3QM r* \ 1/2 _ 3 
Uw ~ sin 2 + 61 VoT r ) ^ ln(2/0 o ) 10- 4 M o /yr I m* 10R o J ° m ' 

(5) 

A test of this model is the presence of very compact H II regions, embedded in the 
outflows of massive protostars. These regions will initially be elongated along the 
outflow axis, but should eventually ionize the entire flow as the star's luminosity 
grows and the accretion and outflow rates finally diminish. This transition 
should be quite rapid for a flow density distribution declining as (rsin#) -2 . 

Of the hot cores in Fig. 1, No centimeter emission is seen (to ~ 0.5 mJy) 
from IRAS 23385+6053 and G34.24+0.13MM, perhaps because of their lower 
luminosities and further distances. The Orion hot core contains a very compact 
thermal source ("I"), discussed below. W3(H20) contains a narrow radio jet 
(Wilner et al. 1999). Shepherd & Kurtz (1999) and Zhang et al. (2002) report 
radio and SiO jets aligned with outflows in G192. 16-3.82 and AFGL 5142. See 
also Garay & Lizano (1999). 
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Figure 2. Radio continuum spectrum of source "I" in the Orion hot 
core. Data at <86 GHz are tabulated in Plambeck et al. (1995). Ad- 
ditional data and limits at 43.1, 98, and 218 GHz are from Chandler 
& Wood (1997), Murata et al. (1992) and Blake et al. (1996), respec- 
tively. Models are described in the text. 

5. The Orion Hot Core and Source "I" 

At a distance of 450 pc, the Orion hot core is the closest example of a massive 
star in formation. The core is self-luminous (-Lbol ~ 1 — 5 x 10 4 L , Gezari et 
al. 1998 [G98]; Kaufman et al. 1998). A weak radio continuum source ("I") 
(e.g. Menten & Reid 1995), located within a few arcseconds of the core center, 
as traced by dust and gas emission (Wright, Plambeck & Wilner 1996), almost 
certainly pinpoints the location of the massive protostar. Note that the Becklin- 
Neugebauer (BN) object (L ~ 10 4 L Q ; 0.02 pc to the NW in projection) is likely 
to be a runaway B star (Plambeck et al. 1995; Plambeck 2002, private comm.), 
that did not form close to the hot core; radio source "L" (0.007 pc to the SW in 
projection; sometimes referred to as "n") is not particularly embedded or lumi- 
nous (G98) and may be an intermediate mass protostar. Note also that bright 
near and mid-infrared features (e.g. IrC2) often result from inhomogeneous ex- 
tinction (G98). X-ray observations (Garmire et al. 2000) can detect lower-mass 
protostars to Ay ~ 60 and reveal several sources (including "L") within about 
0.02 pc projected distance of source "I", which itself is not detected. 

A large scale, wide-angle bipolar outflow extends to the NW and SE of 
the core (Chernin & Wright 1996, and references therein). These authors have 
modeled the flow as being inclined at 65° to our line of sight. At 22 GHz, 
source "I" appears elongated (0".145 by < 0".085, Menten 2002, private comm.) 
parallel to the large scale outflow axis. In a perpendicular direction, SiO emission 
forms a "bow-tie" feature centered on "I" , that may be an inclined or flared disk 
(Wright et al. 1995). SiO should be abundant where dust grains are destroyed. 

We model the thermal radio emission from source "I" assuming the ionized 
gas has a density given by equation (|5|) for two jets of fixed half-opening angle 
9i = 10° and length li = 30 AU, inclined at 65° to our line of sight. These dimen- 
sions are consistent with the 22 GHz size of 0".145 and the lower frequency fluxes 
being due to optically thick emission at temperature Tj = 10 4 K. We investigate 
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models with 6q = 6i and m* = 1,2,4 x lO _4 M0yr _1 . Fig. 2 lists other parame- 
ters. A consistent, but not unique, model accounts for the radio spectrum with 
an outflow from a 2OM protostar accreting at to* ~ 2 x 10~ 4 M Q yr _1 . More 
realistic geometries of the ionized region, derived from the predicted ionizing 
flux, are considered in a future paper (Tan k McKee 2002, in prep.). 

6. Conclusions 

The collisional model for star formation can account qualitatively for some of 
the observational properties of massive stars, such as their tendency to form in 
the centers of clusters. However, it is difficult to achieve the necessary stellar 
densities (~ 10 8 pc -3 ) for this process to be efficient. The standard accretion 
model, modified to account for the high pressures and turbulent, nonthermal 
conditions of massive star-forming clumps can achieve the high-accretion rates 
necessary to overcome radiation pressure and achieve short formation timescales. 
Crowding is not a serious problem. Nearby massive star-forming regions show 
signatures of disks and collimated outflows, suggesting the accretion picture is 
relevant to the formation of stars with masses up to at least ~ 20 — 30 M . 

Acknowledgments. This article summarizes work carried out in collabo- 
ration with Christopher F. McKee. I thank Dick Plambeck and Karl Menten for 
sharing unpublished results on Orion. My research is supported by a Spitzer- 
Cotsen fellowship from Princeton University. 
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